Local quantum critical point in the pseudogap Anderson model: finite-T dynamics 

and uj/T scaling 
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The pseudogap Anderson impurity model is a paradigm for locally critical quantum phase tran- 
sitions. Within the framework of the local moment approach we study its finite-T dynamics, as 
embodied in the single-particle spectrum, in the vicinity of the symmetric quantum critical point 
(QCP) separating generalized Fermi-liquid (Kondo screened) and local moment phases. The scaling 
spectra in both phases, and at the QCP itself, are obtained analytically. A key result is that pure 
oj/T-scaling obtains at the QCP, where the Kondo resonance has just collapsed. The connection 
between the scaling spectra in either phase and that at the QCP is explored in detail. 

PACS numbers: 75.20.Hr, 71.10.Hf, 71.27.-|-a 



I. INTRODUCTION 

Proximity to a magnetic quantum critical point is 
one mechanism for non-Fermi liquid behaviour in heavy 
fermion metals, a topic of great current interesti^^-. Ex- 
perimentally, the temperature below which certain heavy 
fermion materials become antiferromagnetically ordered 
can be tuned to T = either chemically by doping (e.g. 
CeCug-xM^ with M = Au or AgMi^) or through ap- 
plied pressure (e.g. Celng^) or a magnetic field (e.g. 
YbRh2Si2i). The transition point then becomes a quan- 
tum critical point and leads to non-Fermi liquid be- 
haviour extending over large regions of the phase dia- 
gram. 

An example that has received much attention is 
CeCug-aAu x^i^i^ . For a; < — 0.1 it is a paramag- 
net with properties consistent with a highly renormal- 
ized Landau Fermi liquid. The localized magnetic mo- 
ments provided by the Ce atoms, manifest as a Curie- 
Weiss law at high temperatures, are screened by conduc- 
tion electrons at low temperatures in a lattice analogue 
of the familiar Kondo effect for magnetic impurities in 
metals. For x > hy contrast, coupling between local 
moments wins out and the ground state has antiferro- 
magnetic long-ranged order. 

The observed non-Fermi liquid behaviour of 
CeCus.gAuo.A^, probed by neutron scattering ex- 
periments, cannot however be understood in terms 
of the standard Hertz-Millis mean-field description of 
a magnetic quantum phase transition based on long 
wavelength fluctuations of the order parameteriiSiiiii^ 
— the spin dynamics are found to display fractional 
exponents throughout the Brillouin zone, and satisfy 
pure w/T-scaling. This suggests the emergence at the 
QCP of local moments that are temporally critical, 
and that temperature is the only energy scale in the 
problem, cutting of critical quantum fluctuations after a 
correlation time on the order of h/k^T. 

An alternative picture, where critical local fluctuations 
coexist with spatially extended ones, has recently been 
proposed to explain the above results. Si et al^a have 



shown that such a local quantum critical point arises in 
the Kondo lattice model through competition between 
the Kondo exchange coupling — describing the interac- 
tion of conduction electron spin and local moment at 
each lattice site — and the RKKY interaction, which 
characterizes the coupling between local moments at dif- 
ferent lattice sites. At this local QCP the lattice system 
reaches its magnetic ordering transition at precisely the 
same point that the Kondo screening becomes critical. 
Thus the Kondo lattice model, and the self-consistent 
Bose-Fermi Kondo model to which it reduces with an ex- 
tended dynamical mean- field theory treatment, are the 
focus of much current theoretical researcbi2*i^4i^*iSiiLi&. 

However, significant progress in understanding local 
QCPs has followed from the study of a simpler prob- 
lem: the pseudogap Anderson [Kondo] impurity model 
(PAIM), which serves as a paradigm for locally critical 
quantum phase transitions, and is the focus of this work. 
Proposed over a decade ago by Withoff and Fradkin— , 
the model represents a generalization of the familiar An- 
derson [Kondo] impurity model (AIM)2£ describing a 
single magnetic impurity embedded in a non-interacting 
metallic host (for reviews see Ref. l2lll2l . and itself en- 
joying a resurgence of interest in the context of quantum 
dots and STM experiments:^^'^'*. 

In the PAIM the fermionic host is described by a 
power-law density of states vanishing at the Fermi level: 
p{uj) cx \ijjY with r > (r = recovers the case of 
a metallic host). Many specific examples of pseudogap 
systems per se have also been proposed, including vari- 
ous zero-gap semiconductors^^, one-dimensional interact- 
ing systems2&, and the Zn- and Li-doped cupratesSS*^^. 
And the possibility that pseudogap physics might be 
realized in a quantum dot system has also recently 
been suggested^S. The essential physics of the PAIM 
is by now well understood via a number of tech- 
niques, e.g. perturbative scalingiSi^Si^, numerical renor- 
malization group22i22i2iMi2S4^, perturbative renormal- 
ization group^i, local moment approach2^*^2i224£, large- 
N methodsiSiSLiii^, non-crossing approximation^ and 
bare perturbation theory in the interaction strength U—. 
In contrast to the regular metallic AIM, which exhibits 
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Fermi liquid physics and a Kondo effect for all U [or J in 
its Kondo model limit] , the PAIM displays in general crit- 
ical local moment fluctuations and a destruction of the 
Kondo effect at a finite critical Uc [or Jc]. The quantum 
critical point, with spin correlations that are critical in 
time yet spatially local, separates a Fermi liquid Kondo- 
screened phase {U < Uc) from a local moment phase with 
the characteristics of a free spin-1/2. The key features of 
the model are summarized in §2. 

In an effort to understand better the physics of critical 
local moment fluctuations, the most recent work on pseu- 
dogap impurity models has naturally focused on the QCP 
j-j-ggjg7j36j37j39^ which has been established as a non- Fermi 
liquid interacting fixed point for r < 1 (the upper crit- 
ical dimension of the problem). For example, Ingersent 
and Si^^ have recently shown that for r < 1 the dynami- 
cal spin susceptibility at the critical point exhibits LojT- 
scaling with a fractional exponent, which features closely 
parallel those of the local QCP of the Kondo lattice. A 
very recent perturbative renormalization group study'^'' 
has developed critical theories for the transition, confirm- 
ing that, for < r < 1 the transition is described by an 
interacting field theory with universal local moment fluc- 
tuations and hyperscaling, including oj/T-scaling in the 
dynamics. 

In this paper, using the local moment approach 
(LMA), we develop an analytical description of finite- 
T dynamics in the vicinity of the QCP, as embodied in 
the local single-particle spectrum D(uj). Developed origi- 
nally to describe metallic AIMs at T = Q^^J^^ and sub- 
sequently extended to handle finite temperatures^ and 
magnetic fields^^*^, the LMA has also been applied to 
the PAIM2S*^2iia leading to a number of new predictions 
that have since been confirmed by NRG calculations^^. 
We also add that the approach can encompass lattice- 
based models within dynamical mean-field theory, e.g. 
the periodic Anderson model appropriate to the param- 
agnetic metallic phase of heavy fermion materialsi^ii^SiS. 

The paper is organized as follows. After the necessary 
background, section 2, we present in section 3 a number 
of exact results for the scaling spectra in both phases 
of the model, and at the QCP itself, using general scal- 
ing arguments. In section 4 a brief description is given 
of the finite-T local moment approach to the problem, 
along with a numerical demonstration of scaling of the 
single-particle dynamics. An analytical description, aris- 
ing within the LMA, of the finite-T scaling spectrum for 
both phases and at the QCP is given in section 5, and 
is found to be in excellent agreement with the numer- 
ics. The connection between the scaling spectra in either 
phase and that at the QCP itself is explored fully, uj/T- 
scaling of D{u!) at the QCP being one key result. The 
paper concludes with a brief summary. 



II. BACKGROUND 

The Hamiltonian for an AIM is given in standard no- 
tation by 

(2.1) 

where the first term describes the non-interacting host 
with dispersion and density of states p{uj) — J2k ^{^~ 
Efc) (with oj — the Fermi level). The second term de- 
scribes the correlated impurity, with energy Cj and local 
interaction U, and the third the host-impurity coupling 
(via Vik = V). 

For a conventional metallic host the Fermi level den- 
sity of states p(0) 0, and low-energy states are al- 
ways available for screening. In consequence the impu- 
rity spin is quenched and the system a Fermi liquid for all 
U (see e.g.21). For the pseudogap AIM (PAIM) by con- 
trast the host density of states is soft at the Fermi level, 
p{uj) oc (r > 0}iS. Due to the depletion of host states 
around the Fermi level, the PAIM exhibits a quantum 

phase transitionl9i27,28,30,31,32,33,34,35,36,38,39,40^41,42 ^ 

critical U = Uc{r), the quantum critical point (QCP) 
separating a degenerate local moment (LM) phase for 
U > Uc from a 'strong coupling' or generalized Fermi liq- 
uid (GFL) phase, U < Uc, in which the impurity spin is 
locally quenched and a Kondo effect manifest. As r ^ 
the critical Uc{r) diverges (- i/^)33,34,35,38,39,40^ symp- 
tomatic of the absence of a transition for the conventional 
metallic AIM corresponding to r = 0, and whose be- 
haviour is recovered smoothly in the r ^ Iimi1i22i22iiii. 
In the present paper we focus on the so-called symmet- 
ric QCP that arises in both the particle-hole (p-h) sym- 
metric and asymmetric PAIM^i. To that end we con- 
sider explicitly the p-h symmetric model, with ej = —^U 
(and impurity charge ni — '^^{nia) = !)• In this case 
the QCP is known to separate GFL/LM phases for all 
Q ^ J. < 133,34,35,38,39 (foj. ^ > 1 the LM phase alone 

arises for any [/ > 0). 

Our focus here is on single-particle dynamics at finite- 
T, embodied in the (retarded) impurity Green func- 
tion G(w,T) (^ -ie{t){{c;a{t),4„})), with D{uj,T) - 
— ilmG'(ijj, T) the single-particle spectrum. G{lo, T) may 
be expressed as 

G'(tJ,r) [tj+ - A(tj) - I](w,T)]-i (2.2) 

with UJ+ =uj + iO+ and S](w, T) = E^(tj, T) - iY}{uj, T) 
the conventional single self-energy (defined to exclude the 
trivial Hartree term which precisely cancels e-i — —^U), 
such that Y.{uj,T) = -[E(-a;,T)]* by p-h symmetry. 
All effects of host-impurity coupling at the one-electron 
level are embodied in the hybridization function A(a;) = 

EfcKfc['-'+ - efc]"' = ^r(^) - iAi(w)- For the PAIM 
considered here, lS.i{uj) — t:V'^p{ljj) is given explicitly by 
Ai(w) = Ao(|w|/An)'^6'(D - with Aq the hybridiza- 
tion strength, D the bandwidth and 6{x) the unit step 
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function. Throughout the paper we take Aq = 1 as the 
energy unit, i.e. 



Aiito) = \L0\-d{D~\L0\). 



(2.3) 



The real part of the hybridization function fohows simply 
from Hilbert transformation, 



Ar(w) = -sgn(w) 



/3(r)AiH + 



(2.4) 



with /3(r) — tan(^). That this form for the hybridiza- 
tion is simplified is of course irrelevant to the low-energy 
scaling behaviour of the problem, in the same way that 
the usual 'flat band' caricature of the metallic host is im- 
material to the intrinsic Kondo physics of the metallic 
AIM21. 

We now summarize key characteristics of the problem 

at T = 013,27,28,30,31,32,^3,34,35,36,38,39,40,41,42 ^^^^ 

quired in the remainder of the paper, including the essen- 
tial manner in which the underlying transition is directly 
apparent in single-particle dynamics. 



A. GFL phase 



For all U < Uc{r) in the GFL phase (0 < r < i), the 
leading low-o; behaviour of D(lu,0) is entirely unrenor- 
malized from the non-interacting limit and given by 



7rD(tj,0) 



2 ^ " ' M I- 

cos-^(y)|w| 



(2.5) 



This is an exact result^. It arises because T.^/\lu,0) 
vanish as w — > more rapidly than the hydridization 
(oc Iwj''), reflecting in physical terms the perturbative 
continuity to the non-interacting limit that in essence 
defines a Fermi liquid. The leading low-w behaviour of 
S^(ti',0) is in fact given by 



£"(c^,0) 



1 ui 



(2.6) 



as one would cursorily expect from p-h symmetry 
(5]^(w,0) — — 1]^(— tJ, 0)); with the quasiparticle weight 
(or mass renormalization) Z defined as usual by Z = 



[l-(5I]«(c^,0)/ac^), 



By virtue of equation (2.5) 



the most revealing expose of GFL dynamics at T = 
lies in the modified spectral function !F{uj, o)Si^2i22ii&, 
defined generally by 



T{uj,T) = TT sec'^ {^)\uj\'-D{uj,T) 



(2.7) 



such that J-{0, 0) = 1 (i.e. the T — modified spectrum 
is 'pinned' at the Fermi level, a result that reduces to 
the trivial dictates of the Fricdel sum rule for the r = 
metaUic AIM^). 

The Kondo resonance is directly apparent in 

J^(t^, 0)35,38,39,40 _ g_g^ gg^j.g p^gf_ |3| g^jj^ 

figure 3 below. Indeed on visual inspection it is barely 



distinguishable from its r = metallic counterpart, to 
which it reduces for r — 0. The resonance is naturally 
characterised by a low-energy Kondo scale wk, usually 
defined in practice by the width of the resonance^^. Most 
importantly, as ?7 ^ [4— and the GFL— ^LM transition 
is approached, the Kondo scale becomes arbitrarily small 
(vanishing at the QCP itself, where the Kondo resonance 
collapses). For r ^ 1 it is well-established from previ- 
ous LMA work and NRG calculationgSiSi^aia that ljk 
vanishes according to wk oc (1 — U/Uc{r))~ . 

In consequence, in the vicinity of the transition, 
^^(^^,0) exhibits universal scaling behaviour in terms of 
uj/luk^^^- It is of course in this scaling behaviour 
that the underlying quantum phase transition is directly 
manifest^^iili. 



B. LM phase 

The local moment phase is more subtle than naive ex- 
pectation might suggest. Here it is known, originally 
from NRG calculations^^, that the leading low-w be- 
haviour oi D{llj,0) is 



(2.8) 



(with ReG(w,0) sg-n{uj)Tr f3{r)D{uj,0) following di- 
rectly from Hilbert transformation). From this alone, 
simply by inverting equation (2.2), the leading low-cj be- 
haviour of the single self-energy 11^ {uj, 0) is given by 



S^(tJ,0) '"i^" sgn(u;)sin(7rr) 



2ttD{uj,0) 



oc It^r'- (2.9) 



(and likewise 'S}{uj,0) oc S^(w,0)). This divergent low-w 
behaviour is radically different from the simple analytic- 
ity of its counterpart equation (2.6) in the GFL phase, 
and illustrates the basic difficulty: the need to capture 
within a common framework such distinct behaviour in 
the GFL and LM phases, indicative of the underlying 
phase transition and as such requiring an inherently non- 
perturbative description. We do not know of any theoret- 
ical approach based on the conventional single self-energy 
that can handle this issue. 

There is however nothing sacrosanct in direct use of 
the single self-energy, which is merely defined by the 
Dyson equation implicit in equation (2.2). Indeed for 
the doubly-degenerate LM phase, a physically far more 
natural approach would be in terms of a two-self-energy 
description (as immediately obvious by consideration e.g. 
of the atomic limit — the 'extreme' LM case). Here the 
rotationally invariant G{(jJ,T) is expressed as 

G{u;,T) = ^[G^{oj,T) + Gi{u;,T)] (2.10) 

with the Go-(ii',r) given by 

G^iiv, T) = [w+ - A(c^) - t^Lo, r)]-i (2.11) 



4 



in terms of spin-dependent self-energies T,cr{i^,T) (= 
E^(t^,r)-iE^(a;,r)); satisfying 

E,(w,T) = -[E_,(-^,T)]* (2.12) 

for the p-h symmetric case considered. It is pre- 
cisely this two-self-energy framework that underlies the 
LMA3Si39i40,45,46,47,48,49,50 (^j^j^ ^he single self-energy 

obtained if desired as a byproduct, following simply from 
direct comparison of equations (2.10, 11) with equation 
(2.2)); requisite details will be given in section 4. The 
LM phase is then characterised by non-vanishing, spin- 
dependent 'renormalized levels' E^(0,0) = -S5^(0, 0) 
(so called because they correspond to = Ci + S^ (0, 0) 
if the Hartree contribution is explicitly retained in the 
T,^{lo,0)). From this, using equations (2.10,11), the 
low-w behaviour embodied in equation (2.8) (and hence 
equation (2.9) for the single self-energy) follows directly, 
namely 

7rD{cu,0) '"^°|c.|7[S«(0,0)]2 (2.13) 

(assuming merely that ^l^i^i 0) vanishes no more rapidly 
than the hybridization cx; \uj\^). As C/ — > Uc+ and 
the LM^GFL transition is approached, the renormal- 
ized level E^(0,0) ^ (&^; and 2^(0,0) = re- 
mains throughout the ('symmetry unbroken') GPL phase 
U < ^^^38,39.40^ such that the characteristic low-w spectral 
behaviour equation (2.5) is likewise correctly recovered. 

The underlying two-self-energy description is thus ca- 
pable of handling both phases simultaneously, and hence 
the transition between them. Further, since the renor- 
malized level |S^(0, 0)1 vanishes as the LM^GFL tran- 
sition is approached, it naturally generates a low-energy 
scale wl characteristic of the LM phase (defined precisely 
in section 5.2); in terms of which LM phase dynam- 
ics in the vicinity of the transition have recently been 
shown^SiiSi to exhibit universal scaling. This behaviour 
is of course the counterpart of the Kondo scale wk and 
associated spectral scaling in the GFL phase, but now 
for the LM phase on the other side of the underlying 
quantum phase transition. 

III. SCALING: GENERAL CONSIDERATIONS 

Close to and on either side of the quantum phase tran- 
sition, the problem is thus characterised by a single low- 
energy scale, denoted generically by (the Kondo scale 
for the GFL phase; for the LM phase). Given this, and 
without recourse to any particular theory, we now show 
that general scaling arguments may be used to obtain 
a number of exact results for the finite-T scaling spec- 
tra, in both the GFL and LM phases and at the QCP 
itself (where = identically). These are important 
in themselves, and in addition provide stringent dictates 
that should be satisfied by any approximate theoretical 
approach to the problem. 



We consider first the case of T = 0, for which gen- 
eral scaling arguments made by us ini^ serve as a start- 
ing point. In the following we denote the spectrum by 
D{U ; w), with the [/-dependence temporarily explicit. As 
the transition is approached, u = |1 — U/Uc{r)\ — > 0, the 
low-energy scale vanishes, as 

(3.1) 

with exponent a. D{U ; lj) can now be expressed generally 
in the scaling form ttD{U;u!) — u^"-^'^ /u"") in terms 
of two exponents a and b (and with a — GFL or LM 
denoting the appropriate phase), i.e. as 

7rwjL>([/;w) = *„(tj/w,) (3.2) 

with the exponent a eliminated and the w-dependence 
encoded solely in tj/oj* = a). Equation (3.2) simply 
embodies the universal scaling behaviour of the T = 
single-particle spectrum close to the QPT. Using it, three 
results follow: 

(i) That the exponent 6 = r for the GFL phase (and for 
all r where the GFL phase arises). This follows directly 
from equation (3.2) using the fact that the leading uo ^ Q 
behaviour of D{U ; u) throughout the GFL phase is given 
exactly by equation (2.5), together with the fact that 
^Q,(x) is universal; i.e. 

TTUjlD{U;uo)^^^{uo/u,). (3.3) 

(ii) Now consider the approach to the QCP, cj* ^ 
and hence w/w* — > oo. Since the QCP itself must be 
'scale free' (independent of w*), equation (3.3) implies 

vI/„(c.M) 1"!/-^^°° C(r)(|c.|M)-'- (3.4) 

with C(r) a constant (naturally r-dependent). From 
equations (3.3,4) it follows directly that precisely at the 
QCP (w* = 0), 

nD{U,-u) = C{r)\Lo\'"- (3.5) 

which gives explicitly the w-dependence of the {T = 0) 
QCP spectrum. 

(iii) Assuming naturally that the QCP behaviour is 
independent of the phase from which it is approached, 
equation (3.4) applies also to the LM phase. From equa- 
tions (3.2) and (3.4) it follows that the exponent 6 = r for 
the LM phase as well (as one might expect physically). 

The above behaviour is indeed as found in practice 
from the LMA at T = o3**-39-40^ ^-^^^^ ^-^^ scahng for equa- 
tion (3.3) also confirmed by NRG calculations'^^. But the 
important point here is that these results are general, in- 
dependent of approximations (be they analytical or nu- 
merical), and as such holding across the entire r-range for 
which the symmetric QCP exists, i.e. < r < We 
add morover that while the argument above for b = r 
leading to equation (3.3) is particular to the symmet- 
ric PAIM considered explicitly here, recent LMA results 
for the asymmetric PAIM for < r < l'"^, and NRG 
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resultsS^ for the asymmetric pseudogap Kondo model for 
i < r < 1, are also consistent with the exponent b — r. 
The QCP behaviour equation (3.5) then follows directly, 
and the LM A/NRG results of Ref . I2I40I are hkewise con- 
sistent with it. 

The arguments above^S can now be extended to finite 
temperature. The scaling spectra may again be cast gen- 
erally in the form nD{U]iv,T) = T/m°), 
i.e. as 



(3.6) 



expressing the fact that ujID{U ; lu, T) is a universal func- 
tion of cl; = (jj/u)i, and T = T/w*. From this follow three 
results, new to our knowledge:- 

(i) Consider first the approach to the QCP at finitc- 
T, i.e. u) 00 and T — > 00 such that uj/T = uj/T 
is fixed. Again, since the QCP must be scale free, 
^^(w oo,f 00) must be of form fPS{uj/T) with 
the exponent p ~ —r from equation (3.6), i.e. 



(3.7) 



(ii) From equations (3.6,7) it follows directly that pre- 
cisely at the QCP (w* = 0): 



TrT^D{U,;io,T)^ S{uj/T) 



(3.8) 



i.e. T'^ D[Uc\ijJ,T) exhibits pure /T -scaling. 

(iii) The large- a; behaviour of S{x) — and hence the 
'tail' behaviour of the QCP scaling spectrum S{uj/T) — 
may be deduced on the natural assumption that the lim- 
its cj* and T ^ commute; i.e. that the T = Q QCP 
spectrum may either be obtained at T = from the limit 
— > (as in equation (3.5)), or from the T — > limit 
of the finite-T QCP scaling spectrum equation (3.8) (in 
which — > has been taken first). With this, compari- 
son of equations (3.5,8) yields the desired result 



\/T- 



Cir)i\co\/T)- 



(3.9) 



The behaviour embodied in equations (3.7-9) follows 
on general grounds. In particular, the w/T-scaling of 
the QCP scaling spectrum is a key conclusion. Nei- 
ther is such scaling specific to single-particle dynamics. 
For the dynamical local magnetic susceptibility, recent 
NRG calculations'^^ are consistent with w/T-scaling at 
the QCP (for < r < 1), as supported further by ap- 
proximate analytic results for small r based on a proce- 
dure analogous to the standard e-expansion. The preced- 
ing arguments do not of course determine the functional 
form of the scaling spectra embodied in the ^>a{oJ,T) 
and S{uj/T), save for their asymptotic form equations 
(3.4,9): for that a 'real' (inevitably approximate) theory 
is required, as considered in the following sections; and 
from which the general results above should of course 
be recovered. We note in passing that the arguments 
above also encompass the role of a local magnetic field, 
h — ^gpBH\oc (with h = /i/w* in the following): re- 
placing r by in equations (3.6-9) gives the appropriate 



scaling behaviour in terms of Cj and h (for T = 0). At 
the QCP in particular D{Uc]OJ,h) exhibits pure uj/h- 
scaling, 



Trh'^D{Uc;uj,h) = S'{uj/h) 



(3.10) 



with the large- 1 w I //i 'tail' behaviour S'{uj/h) 
C{r)i\u;\/h)-\ 



A. The low-energy scale 

Before proceeding we return briefly to an obvious is- 
sue: the physical nature of the low-energy scale w*. In 
the GFL phase for example, how is it related to the quasi- 
particle weight Z = [l — {d'E^{uj; 0)/dLLj)ui=o]~'^ , which we 
likewise expect to vanish as U Uc{r) — , 



Z 



(3.11) 



with exponent ai? That question may be answered most 
simply (albeit partially here) by considering the limit- 
ing low-frequency 'quasiparticle form' for the impurity 
Green function, whose behaviour reflects the adiabatic 
continuity to the non-interacting limit that is intrinsic to 
the GFL phase. This follows from the low-w expansion 
of the single self-energy, S(w;0) ~ — (^ — l)uj (with the 
imaginary part neglected as usualSi); using which with 
equations (2.2-4) yields the quasiparticle behaviour 



TrD{U;uj) 



(|+sgn(c.)/3(r)|a;|'-)2 + |c.|2'- 



(3.12) 



which itself must be expressible in the scaling form equa- 
tion (3.3). This in turn implies that and Z are related 

by 



= [kiZ] i--- 



(3.13) 



(with ki an arbitrary constant up to which low-energy 
scales are deflned), such that equation (3.12) reduces to 
the required form 



7Tcu:DiU;io) 



sgn{Lu)(3{r)\Cjn^ + \Cj\^^ 



(fcipl 

(3.14) 

with Lu — ijjjLo^. For the GFL phase the uj^ (= lok) 
scale and the quasiparticle weight Z are thus related by 
equation (3.13) (and hence the exponents a (equation 
(3.1)) and a\ (equation (3.11)) by a\ — a(l — r)). That 
behaviour is correctly recovered by the LMA, which in 
addition shows to be equivalently the characteristic 
Kondo spin-flip scale associated with transverse spin ex- 
citations; see section 5.1 below. 

In contrast to the GFL phase where the leading low- 
Lo spectral behaviour is D{U;uj) oc \uj\~^ as above, the 
corresponding behaviour for the LM phase is given by 
ttD{U]lo) cx |w|7[f:^(0;0)]2 (equation (2.13)) in terms 
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FIG. 1: Standard LMA diagram for the dynamical contribu- 
tion Eo- to the self-energies. 



of the 'renormalized level' T,^ {0; 0), which is non-zero in 



the LM phase and vanishes as U ^ Uc{r)- 



I.e. 



(3.15) 



with exponent 02. Since equation (2.13) must be ex- 
pressible in the scaling form equation (3.3) it follows that 
uj^,{= lul) in the LM phase is related simply to the renor- 
malized level by 



[fc2|S^(0;l 



(3.16) 



(and hence the exponents a and 02 by 02 — ar). We 
emphasize again that the exponents obtained above are 
exact, valid for all < r < i. 

After the above general considerations we now consider 
the finite-temperature scaling properties of the PAIM 
within the local moment approach. 



IV. LOCAL MOMENT APPROACH: FINITE T 

Regardless of the phase considered the LMA uses 
the two-self energy description embodied in equations 
(2.10,11), with the self-energies separated for convenience 
as 

t„{ij+,T) = -^U\fi\+E^{uj+,T) (4.1) 

into a purely static Fock piece with local moment \fi\ 
(that is retained alone at mean- field level); and an all 
important dynamical contribution 'Ei^{uj~^ ,T) containing 
the spin-flip physics that dominates at low-energies. The 
Ga{uj,T) (equation (2.11)) may be recast equivalently as 

G^iuj, T) - [g^Lo)-^ - T)] (4.2) 

in terms of the mean-field (MF) propagators 

g^(^) = [c.+ + |c/H-A(c.)]-i (4.3) 

and with Yi„ = Ti„[{QcY\ a functional of the {Qc{ij)}. 

The standard LMA diagrammatic approximation for 
the i-ordered two-self-energies is depicted in figure 1. For 
full details, including their physical interpretation, the 



reader is referred to2S4SiH. In retarded form at finite 
temperature, S|(<j-') is given bj*^ 

st(c.+,t) = r ^ r dL.2 x^-{^i.T) 

J-ao J-00 



^1(^2) 



o; + 0^1 — 0^2 + i'7 



3('^i;'^2) 



(4.4) 



where D^iuj) = —-lmQ„{uj) is the MF spectral density; 
and giLUi;L02) = el^Lo^) [1 - /(c^2;T)] + 0(c^i)/(c^2; T) 
with /(w;T) = [1 + exp{uj/T)]'^ the Fermi function. 
X~^~{co,T) = ImH"' (u), T) is the spectrum of transverse 
spin excitations, with H"' (w,r) the finite-T polariza- 
tion propagator, calculated via an RPA-like p-h ladder 
sum/®. Sxl*^) = ~[^t(~'^)]* follows trivially by particle- 
hole symmetry; we thus choose to work with I]|(w). 

The T-dependence of the ^^id also of the 

local moment = \fJ.{T)\ entering the MF propaga- 
tors (equation(4.3)) and self-energies (equation (4.1)), is 
found to be insignificantly small, provided we remain 
uninterested in temperatures on the order of the non- 
universal energy scales in the problem. This is pre- 
cisely as for the r = case^ and we thus work with 
X+"'(w,T) ~ x+"('^,r = 0) = Imn+-(w) and local mo- 
ment |/i(0)| = \fj,\. The remaining T-dependence, de- 
scribing the universal scaling regime, is controlled solely 
by g{uji]uj2) entering equation (4.4). 

In describing the GFL phase within the LMA, the key 
concept is that of symmetry restoration (SR)^242i^^*^: 
self-consistent restoration of the symmetry broken at 
pure MF level. In mathematical terms this is encoded 
simply in S|(0+,T = 0) = S|(0+,T 0) and hence, via 
particle-hole symmetry. 



Et(0+;T = 0) = 



ST(0)-ic/|A^| 



0. 



(4.5) 



Equation (4.5) guarantees in particular the correct Fermi 
liquid form !F{0,0) = 1 (see equation (2.7)). It is satis- 
fied in practice, for any given U < Uc, by varying the 
local moment |/^| entering the MF propagators from the 
pure MF value |/io| < \lA until equation (4.5) — a sin- 
gle condition at the Fermi level w = — is satisfied. In 
so doing a low-energy scale is introduced into the prob- 
lem through a strong resonance in ImH"' {uj) centred on 
Lj — Lj„^ = Lo-^{r). This is the Kondo or spin-flip scale 
LJa\ oc wk and sets the timescale for symmetry restora- 

tion, r ~ /t/t^„i38,40,45,47^ 

In the LM phase, by contrast, it is not possible to sat- 
isfy the symmetry restoration condition equation (4.5) 
and — |/ion'*'^°: Imll^ (cj) then correctly contains 
a delta-function contribution at cj = 0, reflecting phys- 
ically the zero-energy cost for spin-flips in the doubly- 
degenerate LM phase. The 'renormalized levels' (see sec- 
tion 2.2) {So-(0, 0)} are here found to be non-zero and 
sign-definite. 

The GFL/LM phase boundary, i.e. Uc = Uc{r), may 
then be accessed either as the limit of solutions to equa- 
tion (4.5), where lo^ — > 0; or, coming from the LM side. 




-2 -10 12 



FIG. 2: GFL phase scaling spectrum T{uj,T) — 
7rsec^(^r)|ti;|''Z)((jj, T) versus oj — co/ujm, for fixed T — 
T/uJm ~ 1 and four values of the interaction strength U ap- 
proaching Uc — 16.4. The inset shows T{uj,T) on an absolute 
scale. 



as the limiting U for which the {Eo-(0,0)} vanish. In ei- 
ther case the same Uc{r) obtains; and both GFL and LM 
phases are correctly found to arise for all < r < ^ , while 
solely LM states occur for all r > i and U > 0^. For 
r — > the LMA recovers the exact result that Uc = S/irr 
which, together with the low-r result a — l/r (see equa- 
tion (3.1))"^^ recovers the exact exponential dependence 
of the regular r = metallic AIMili. 

The practical strategy for solving the problem at finite- 
T is likewise straightforward. Once the local moment 
is known from symmetry restoration (for any given 
r and U), the dynamical self-energy S]|(ijj+,T) follows 
from equation (4.4). The full self-energy then follows 
immediately from equation (4.1) and G{lli,T) in turn via 
equations (4.2) and (2.10). We now give numerical results 
obtained in this way for the GFL phase at finite-T. 



A. Numerical results: GFL phase scaling 

Here we show that the LMA correctly captures scal- 
ing of the single-particle spectrum D{uj,T) as the tran- 
sition is approached, U — > Uc{r) — . That is, ljI-^^D{uj,T) 
scales in terms of both a) — w/wm and T — T/wm, as ex- 
pressed in equation (3.6); or equivalently that T{lo,T) = 

Figure 2 shows the numerically determined spectral 
function J-{uj,T) versus Co = w/w,n for r = 0.2, a fixed 
temperature T — T /lj^ — 1 and four values oiU < Uc{r) 
as the transition is approached {Uc — 16.4). The spectra 
— very different on an absolute scale as shown in the 
inset to the figure — indeed clearly collapse to a common 
scaling form as C/ — > Uc{r)— and Wm becomes arbitrarily 
small. 

This behaviour is not of course particular to the choice 
T = 1, but rather is found to arise for all (finite) T. 
Figure 3 shows the resultant scaling spectra for a repre- 




()i 1 < ' < 1 < ' < ' 1 

^ -40 -20 20 40 

UJ 

FIG. 3: Numerical results: GFL phase scaling spectrum 
J-{ijj,T) versus uj = oj/um for a range of T = T/ujm- The 
generalized Kondo resonance is thermally destroyed, with the 
principal effects of temperature arising for \Cb\ < T. The inset 
shows uj^Diuj, T) versus a) on a logarithmic scale for the same 
T, increasing top to bottom; see text for comments. 



sentative range of T, labelled in the figure. For T = 
the GFL phase Kondo resonance is fully present and 
J^{uj = 0,T = 0) = 1, as discussed in section 3. As 
temperature is raised through a range on the order of 
the Kondo scale itself (T ~ 0(1)) the Kondo resonance is 
thermally destroyed. First it is split — reflecting the fact 
that the Fermi level D{iLj = 0,T) is finite for T > (see 
inset) — and then progressively eroded as T increases. 
For any given temperature T the principal effect of tem- 
perature arises on frequency scales \Lij\ < T, such that for 
UJ ^ T the spectrum is essentially coincident with the 
T — limit. This behaviour is observed clearly in the 
inset to figure 3 wherein e.g. the T = 1 scaling spectrum 
is coincident with that for T = for \ui\ > 10; and that 
for f 10 coincides with f = for > 50. We also 
add that the small spectral feature at (D ~ 1, seen e.g. in 
the inset to figure 3, is entirely an artefact of the specific 
RPA form for Imll^ (w) that, as discussed in detail in 
previous work^, can be removed entirely with both little 
effect on the appearance of the spectrum and no effect 
on any asymptotic results. 



V. RESULTS 

We now turn to an analytical description of the finite-T 
scaling spectrum in the strong coupling (= strongly cor- 
related) limit oi U ^ I. As in Ref . l3il and in contrast 
to the results of section 3, this formally restricts discus- 
sion to r ^ 1, although as demonstrated below such an 
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analysis does rather well in accounting for the numerical 
results given in section 4.1 for r = 0.2. For this reason 
we present only analytical results in the remainder of the 
paper, where we focus in turn on the GFL phase, LM 
phase and the QCP itself. 



A. GFL phase 

To obtain the scaling behaviour of D{ijj) we consider 
finite Cj = Lo/um in the formal limit Wm 0, thereby 
projecting out irrelevant non-universal spectral features 
(e.g. the Hubbard bands at cj ^ ±J7/2). Using equations 
(2.10,11) the scaling spectrum w^_D(w,T) is thus given 
by 



H{y) is defined as 



exp(a;) + 1 



V 



1 



1 



uj + y uj~y 



(5.6) 



and is trivial to evaluate numerically; its ?; >> 1 and 
y « 1 asymptotic behaviour is also obtainable in closed 
form, as used below. From equations (5.4,5), we thus 
obtain 



«)-^I]f(c.,T) 



7(r)[|l+^r-l] 
\l + u:\ 



4 ~ 

TT 



T 



(5.7) 



<i?(a,,r) = - 



where (from equations (2.3,4)) the hybridization 
A(a;) = A((I)a;,„) reduces simply to (wj^)^-^ A(ti;) ~ 
— sgn(a;)[/3(r)-|-i]|(i|'' (and, trivially, the bare uj = lo,^uj of 
equation (2.11) may be neglected). We consider explic- 
itly CJ > in the following, since D{—uj,T) = D{uj,T) 
from particle-hole symmetry. 

Our task now is to demonstrate that (tJ^-,)~^S|(a;, T) 
is a function solely of Q and T in the GFL phase. It 
follows from equation (4.4) that 



The complete GFL scaling spectrum now follows from 
equations (5.1) and (5.3,7); and, as required, lu'^D{uj,T) 
is seen to be a universal function of d) = w /uj^ and T = 

Before proceeding to a comparison with figure 3 how- 
ever, we add that uj^D(u},T) is equivalently a universal 
function of cu/T = cu' and T^^. To see this note that 
equations (5.3,7) may be rewritten as 





1 


r 






f 


f 






T 


VT J 



(5.8) 



and 



/oo 
du;ix^-{i^i)Dl{cJi+u;) (5.2) (T^-^Ef (cj, T) 



(with the temperature dependence of /(cj, T) temporarily 
omitted). In the requisite strong coupling regime C/ ^ 1, 
(w,T) ~ Imll"' (w) = TrS{uj — Wm) is readily shown 
to reduce asymptotically to a (5-function centred on lo = 
uj,n; and thus ^{co,T) = nU^D°^{uJn,[l + oj])f{l + Cj^f) 
(noting that f{uj,T) depends solely on the ratio io/T). 
Using equation (4.3), -Dj(cj) ^ 4|a;|'"/7rC/^ in strong cou- 
pling; whence 



-H 





r 




rj^ — r 


UJ 










) 


T 





(5.9) 



(with notation f{z) = f{z, 1) — [exp(z) + 1] ^ for the 
Fermi functions), which depend solely on u' and T^^. 
From equation (5.1) it follows directly that 



«)-^I]|(c., T) = 4|1 + c^r/(l + ^, T) 



(5.3) 



TTT^Diuj^T) = --Im^ {(r'-)-i [a(c.) - S.(tJ,T)]} 

(5.10) 



which is thus of form 



which is indeed universally dependent solely on Cj and T. 
We now consider the real part Yi^ {u), T). For T = it is 
known to be given by^ 



^T-D{uj,T) = Sgfi. 



UJ 1 



(5.11) 



«)~isf(cj,0) = -7(r)[|l+cir-l] 



(5.4) 



where 7(r) — 4/sin(7rr) ~ 4/7rr, scaling solely in 
terms of uj and correctly satisfying symmetry restoration 
(S^(0, 0) = 0) for the GFL phase. The finite-T difference 

5i:f{uj,T) = Sf(w,T) - Ef(a;,0) is readily calculated 
via Hilbert transforms using equation (5.2). After some 
manipulation it is found to be given by 



«)-i<5Ef(c.,T) = -if'-ij('^ 



(5.5) 



The importance of this result is that from it the finitc-T 
spectrum at the QCP itself {cf equation (3.8)) may be 
obtained explicitly, since the QCP corresponds to uj^^ = 
i.e. T^^ = 0. We consider this issue explicitly in section 
5.3. 

First we compare directly the analytic results above 
for the GFL phase, with the numerical results given in 
figure 3. Figure 4 shows the modified spectral function 
J-{uj) versus uj = uj/uj^ for r = 0.2 arising from equa- 
tion (5.1,3,7), for the same range of 'scaled' temperature 
T = T/ujai. The agreement is seen to be excellent: all 
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the high-frequency asymptotic 'tails' equation (5.13) is 



7I"<n-C(w,r) 



16 



(5.14) 



precisely as found hitherto for T ~ 0^; and since the 
scale u!^^ drops out of equation (5.14), the low-w be- 
haviour of the T = QCP spectrum itself (where uj„i — 
0) follows immediately as^ 



7rD(w,0) 



2^2 



16 



OJ 



(5.15) 



which result is believed to be asymptotically exact as 
r ^ 0. 



FIG. 4: Analytical results: GFL phase scaling spectrum 
J-(u},T) versus ui = to/uini for r = 0.2 and the sequence of 
temperatures shown. 



features of the numerical spectra are captured by the an- 
alytical form; quantitative differences are small, and arise 
because the equations hold asymptotically as r — > 0. 

Two specific points concerning the large-w behaviour 
of the GFL phase scaling spectrum may now be made. 
First, the observation (section 4.1) that for frequen- 
cies ^ T the scaling spectra are effectively inde- 
pendent of T. For ^ max(l,T) the argument 
of H{\1 + Cj\/f) is large and H{\1 + Cj\/f) ~ may 
be neglected in equation (5.7) {H{y) ~ — vr^/ey^ for 
y » 1). Hence «)-iSf(c^,r) - -l{r)[\uY - 1] 
and (ujl,)-'^t\{uj,T) ~ 4|w|''6l(-[H- w]) from equations 
(5.7,3), and from equation (5.1) it follows that 



1 



2\Co\^ i[/3(r)-t-7(r)(l-p|-'-)]2-f 1 



(5.12) 



7(r)(l 



The key point here is that equation (5.12) is independent 
of T, showing explicitly that for \Cj\ ^ max(l,T) the 
scaling spectrum indeed coincides with its T = limifc^?. 

Second, equation (5.12) shows that the leading large Co 
behaviour {\uj\'^ ^ 1) of the GFL scaling spectrum is 



TTUJ]^D{UJ,T) 



2\Q\^ {[P{r)+j{r)]^ + 1 



W - 7W]^ + 25 



(5.13) 



- i.e. uj^-^D{uj, T) oc \Cj\~^ which, for small r, is very slowly 
varying. The onset of this powerlaw behaviour for |ti;|'' ^ 
1 is readily seen in the inset to figure 3. 

Finally, since 7(r) — 4/sin(7rr) ~ 4/7rr and /3(r) = 
tan(^) ~ ^ for r « 1, the leading low-r behaviour of 



B. LM phase 

While the general form equation (5.10) for D{uj,T) 
naturally applies to both phases, symmetry is not of 
course restored for U > Uc in the LM phase22ii£: the 
renormalized level ^^■(O, 0) is non-zero. But it necessar- 
ily vanishes a,s U ^ Uc+ and the LM— >GFL transition is 
approached, and in consequence determines a low-energy 
scale wl characteristic of the LM phase as discussed in 
§3.1. Given explicitly by 



c.L = [|Sf(0,0)|/7M]^ 



(5.16) 



this is the natural scaling counterpart to in the GFL 
phase, such that the T = LM phase spectrum scales 
universally as a function of (D = ui/uji^^^. 

The required results for E|(w,T) in the LM phase are 
obtained most directly from their counterpart equations 
(5.8,9) in the GFL phase (in which l/f = uj^/T and 
l/uj — uj-a^/uj), by setting = therein — recall that 
the spin-flip scale vanishes identically throughout the LM 
phase, where {^iT) ~ Imn"* [uj) = ttS{uj). Equa- 
tion (5.8) then yields 



{T^)-'mu,T) = A\u'\^f{Lo'). 



(5.17) 



T-^[Ef (w,r) - Sf (0,0)] is likewise obtained by setting 
Will = in the right side of equation (5.9). Using equation 
(5.16) for Sf (0,0) = (0,0)1, this gives 

(r'-)-iEf (c.,r) = ~-f{r)[f-^ + Ico'n - -Hito') (5.18) 

TT 

in which T = T/uJi^ (and lu' = uj/T = Oj/T). 

Equations (5.17,18) and (5.10) generate the LM phase 
scaling spectrum. It is thus seen to be of form 



TrT''Diiu,T) = Slm 



UJ 1 

T' T 



(5.19) 



scaling universally in terms of uj' and T; or, entirely 
equivalently, that w£_D(w, T) scales in terms oi lu = uj/uh 
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and T . The scaling spectra so obtained are shown in fig- 
ure 5 for a range of T. It is again readily shown that 
for \uj\ ^ T the LM phase scaling spectra are effectively 
independent of T and coincide with the T = limit'^^, 
and from which the low-w behaviour of the T — QCP 
spectrum (equation (5.15)) correctly follows by taking 



C. Quantum critical point 

We turn now to the quantum critical point itself. It has 
been shown in §3 that the QCP spectrum must exhibit 
pure uj' = Lu/T scaling, ie that ttT''D{uj,T) = S{lu'). In 
§s 5.1,2 we have also shown that this general behaviour is 
correctly recovered by the local moment approach (which 
is distinctly non-trivial, bearing in mind that the QCP is 
an interacting, non- Fermi liquid fixed point). Our aims 
now are twofold. First to obtain explicitly the scaling 
function S{uj') arising within the LMA for r <C 1; and 
then to show that it is uniformly recovered from GFL/LM 
scaling spectrum as T = T/w* — > oo (we remind the 
reader that = uj^ in the GFL phase and = ojl for 
the LM phase). 

Taking = in equations (5.17,18) (or equations 
(5.8,9)) gives 

{Tn-'mLo,T)^A\u;rf{Lo') (5.20) 



(r'-)-isf (c., T) = -7(r)|c.T - -Hiu;'). (5.21) 

TT 

The QCP scaling spectrum ttT''D{uj,T) = S{uj') now 
follows directly from equation (5.10). 

Figure 6 shows the resultant QCP spectrum for r = 
0.2, as a function of u' = ui/T on a logarithmic scale. 



From this it is clear that two distinct power law be- 
haviours dominate. For \uj'\ <C 1, where H{lu') ^ 
r-i(l - IljT) + c (with c = ln(fexp(-C)) ~ -0.125 
a constant, and C Euler's constant), it follows that 

7rT'^i^(c^,T) ^^Ic^T'- (5.22) 
16 

This behaviour is clearly seen in the figure 6, and in prac- 
tice sets in for |u;'| < 1. 

For \uj'\ » 1 by contrast, where H{uj') ~ may be 
neglected, equations (5.10,20,21) give 

nT^Dico,T)^'^'}^^'^\co'r- (5.23) 

This form is also clearly evident in figure 6, in practice 
setting in above \uj'\ > 1. Such behaviour is also precisely 
the general form deduced in section 3 (equation (3.9)) on 
the assumption that the limits a;* — > and T — > com- 
mute; with the exphcit r-dependence C(r) = 37r^r^/16, 
which we believe to be asymptotically exact as r ^ 0, ob- 
tainable because the result now arises from a microscopic 
many-body theory. 

As discussed in section 3 and above, the GFL/LM scal- 
ing spectra may be expressed as functions of w' = lu/T 
for any given T. One obvious question then arises: how 
do the GFL/LM scaling spectra evolve with temperature 
T to approach their ultimate limit of the QCP spectrum 
as T — > oo? This is illustrated in figure 7 where the 
analytic results of sections 5.1 and 5.2 are replotted ver- 
sus Lu' — Lu/T for a sequence of increasing temperatures 
T = T/uj^. The QCP scaling spectrum is shown as a 
thick dashed line. From figure 7 we see that — for all T — 
the GFL and LM scaling spectra coincide both with each 
other and with the QCP scaling spectrum for :i> 1. 
This is a reflection of two facts. First, as discussed in 
sections 5.1 and 5.2, that the high-frequency behaviour 
of the GFL and LM phase spectra coincide with their 
T — Q limits. And second, that the high-frequency tails 
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FIG. 7: The left (right) panel shows the GFL (LM) phase scaling spectrum T'^D{uj,T) = ^gflI^^', T"^) versus lj' = uj/T for 
f = 10"\ 10°, 10\ 10^, 10^, 10*, lO'^ and oo (in sequence top^bottom for the GFL phase and bottom^top for LM). In either 
case the QCP spectrum T = oo is shown as a thick line; the dotted line shows the T spectrum. 



of the r = GFL/LM scaling spectrum coincide with 
those of the T — QCP spectrum. 

Most importantly, however, figure 7 shows that the 
QCP scaling spectrum is obtained 'uniformly' from the 
GFL/LM phase scaling spectrum upon increasing tem- 
perature T. That is, with progressively increasing T 
the GFL/LM phase scaling spectra ever increasingly co- 
incide with the QCP scaling spectrum (thick line) over 
a larger and larger frequency interval, such that in the 
limit T oo the QCP spectrum for all uj/T is obtained 
smoothly. This is an important result, providing as it 
does a direct connection between the QCP scaling spec- 
trum and scaling spectra in the GFL or LM phases at 
any finite T. 

VI. SUMMARY 

We have considered a local moment approach to the 
pseudogap Anderson impurity model, close to the sym- 
metric quantum critical point where the Kondo resonance 
has just collapsed. Building on previous work^2iiM^, 
we have focused on single-particle dynamics at finite 



temperature, obtaining an analytical description of the 
finite-T scaling behaviour for small r, in both generalized 
Fermi liquid (Kondo screened) and local moment phases. 
A key result obtained on general grounds is that pure 
w/T-scaling obtains at the QCP itself, consistent with 
an interacting fixed point and recent results for the lo- 
cal dynamical susceptibility"^^. We have succeeded both 
in obtaining explicitly the QCP scaling spectrum, and 
in understanding its continuous emergence with increas- 
ing T/w* from the scaling dynamics appropriate to the 
Kondo screened and local moment phases on either side 
of the quantum phase transition. 

Related results have also been obtained for the PAIM 
with finite local magnetic field, h — g^^Hioi^^; i.e. 
pure (jj/Zi-scaling of the single-particle spectrum at the 
QCP. These will be discussed elsewhere. 
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